Evolution of unknotting strategies for knots and braids 
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Abstract. This paper explores the problem of unknotting closed 
braids and classical knots in mathematical knot theory. We apply 
evolutionary computation methods to learn sequences of moves that 
simplify knot diagrams, and show that this can be effective both when 
the evolution is carried out for individual knots and when a generic 
sequence of moves is evolved for a set of knots. 



Isotopy of embeddings descends to certain allowable local moves 
on diagrams which were first studied by Reidemeister |16| and by 
Alexander and Briggs [2|. These Reidemeister moves are depicted 
in Figure |2] Two knots or links are isotopic if and only if their dia- 
grams are connected by a finite sequence of Reidemeister moves and 
continuous deformations of the ambient projection plane. 



1 Introduction 
1.1 Knots and links 
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Knot theory is currently one of the richest and most vibrant areas of 
pure mathematics, having connections not only with other topics in 
algebraic and geometric topology, but also with many other branches 
of mathematics, as well as mathematical physics 1201 and biochem- 
istry fW\. 

A full introduction to the study of knots and links is beyond the 
scope of this article, but a readable introduction may be found in, for 
example, the book by Cromwell |i8J, and a more comprehensive but 
still accessible survey in Rolfsen's classic text 1171 . 





Figure 1. Examples of knots: the trefoil (3i), Conway's knot (lln34) and 
the Kinoshita-Terasaka knot (ll?i42) 

We define a knot to be an isotopy class of embeddings K : ^ 
R^, where 5"^ = {e^" £ C : sC < 27r} denotes the standard 
unit circle; informally, this is a set of placements of a closed circle 
in space. A link is a knot with more than one circular component, 
that is, an (isotopy class of an) embedding L: U . . . U ■— > 
R"^. (Alternatively, a knot may be regarded as a link with a single 
component.) 

We generally represent knots and links with diagrams in the plane: 
projections of the embedded circle(s) where each double intersection 
point is equipped with crossing information, and we disallow cusps, 
tangencies or triple-points. Examples may be seen in Figure [T] Iden- 
tifiers such as 3i and 817 refer to the table in Rolfsen's book 1 17], 
while identifiers of the form lla367 and Il7i34 refer to, respec- 
tively, alternating and non-alternating knots in the census of Hoste, 
Thistlethwaite and Weeks llOi . 

^ Mathematics Institute, University of Warwick, Coventry CV4 7AL, UK. 

Email: Nicholas.Jackson@warwick.ac.uk 
^ School of Computing, University of Kent, Canterbury CT2 7NF, UK. Email: 

C.G.Johnson@kent.ac.uk 





«2, 



Figure 2. Reidemeister moves 

There are a number of different measures of the complexity of a 
given knot or link K, the best known of which is the crossing num- 
ber: the minimal nuinber of crossings over all possible diagrams for 
K. Related to this is the unknotting number u{K): the minimal num- 
ber, over all possible diagrams for a knot, of crossings which must 
be changed in order to obtain a trivial knot. The trefoil in Figure [T] 
has unknotting number ii(3i) = 1: it may be seen to be nontrivially 
knotted (that is, not isotopic to an unknotted circle) but changing any 
single crossing results in a diagram which may be transformed (by 
means of an R2 move followed by an Ri move) into an unknotted 
circle. Less obviously, the other two knots in Figure[T]are also known 
to have unknotting number 1 . 

The unknotting number u{K) is a conceptually simple measure 
of the complexity of a given knot K (broadly speaking, the higher 
the unknotting number, the more knotted the knot in question) but 
one which is often not straightforward to calculate. According to the 
Knotlnfo database 1 12 1, the unknotting number is currently unknown 
for nine of the 165 prime knots with ten crossings (and also for many 
knots of higher crossing number); for these nine knots, the unknot- 
ting number is is known to be either 2 or 3, due at least in part to work 
on Heegaard Floer homology by Oszvath and Szabo 1 15] which ruled 
out unknotting number 1. Recent work by Borodzik and Friedl ^ 
has introduced a new invariant which provides a lower bound of 3 
for the unknotting number of twenty-five otherwise difficult cases of 
knots with up to twelve crossings. 

In this paper, our main goal is not necessarily to find optimal 
bounds on the unknotting numbers of currently unresolved cases, but 
to explore the possibilities afforded by applying evolutionary com- 
puting techniques to pure mathematical problems in group theory 



and geometric topology, and to try to obtain some qualitative under- 
standing of the search landscape for these problems. The unknotting 
problem is relatively straightforward to describe and implement, and 
thus provides a good candidate for a preliminary investigation of this 
type. 

1.2 Braids 

The theory of braids was first seriously investigated by Artin 1 3 1 , and 
again a full treatment is far beyond the scope of this article, so we will 
discuss only those aspects essential for what follows. A more com- 
prehensive discussion may be found in either the book by Hansen 1 9| 
or the classic monograph by Birman jSj. 

We define a geometric braid on n strings to be a system of n dis- 
joint, embedded arcs A = {Ai, . . . , A„} in x [0, 1], such that the 
arc Ai joins the point Pi = {i, 0, 1) to the point = (''"(«), 0, 0), 
where r denotes some permutation of the numbers {!,..., n}, such 
that each arc Ai intersects the intermediate plane x {z} exactly 
once, for all < z < 1. Figure [3] shows an example of a 4-string 
braid. 




Figure 3. A 4-string braid 

The elementary braid ai, for 1 ^ i n—1, is the n-string braid 
in which the (i-l-l)st string crosses over the ith string, and no other 
interactions take place; its inverse is the braid in which the ith 
string crosses over the (i+l)st string (see Figure |4]l. Any n-string 
braid /3 may be represented (although not, in general, uniquely) as a 
concatenated sequence of elementary braids. 
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Figure 4. The elementary braids <Ti and 

We consider two n-braids /3i and P2 to be equivalent if they are 
related by an isotopy which keeps their endpoints fixed. In the lan- 
guage of elementary braids, this translates to the following identities: 



(Jiai+iai 



HO"i(Ti + l 



for|i- j| > 1 (1) 
for 1 sC i ^ n - 1 (2) 



Geometrically, the first of these corresponds to moving two non- 
interacting elementary braids past each other, and the second is es- 
sentially the third Reidemeister move Rs- 

We may define the braid group Bn by the following presentation: 



CTl, . . . , cr„_i 1 ^ ■ ^ 1 / 

(yiOi + lCJi = (Ji + -l(JiGi + \ i ^ J ^ n — i 



It can be shown that the group defined by this presentation is isomor- 
phic to the group obtained by imposing the obvious concatenation 



operation on the (in general, infinite) set of all n-string geometric 
braids. In this latter group, the identity element is the trivial n-braid 
(the one with no crossings) and for any braid the inverse /3~^ may 



be obtained by reflecting P in the horizontal plane ] 



There are other, equivalent constructions of the n-string braid 
group, including one in terms of the fundamental group of a particu- 
lar configuration space factored by an action of the symmetric group 
Sn but these will not concern us here. 

Given a braid j3 £ _B„, we can obtain a link /3 by the closure oper- 
ation depicted in Figure[5] that is, we join each point Pi — {i, 0, 1) 
to the point Qi — (i,0, 0). In fact, Alexander's Theorem 11] (see 
also Birman |i5j Theorem 2.1]) states that any knot or link can be 
obtained in this way; an explicit algorithm may be found in the pa- 
per by Vogel f 19]. Note that a closed-braid presentation need not be 
minimal with respect to the crossing number of the knot. That is, an 
71-crossing knot might not have a closed-braid presentation with n 
crossings. Table [T] lists several examples of non-minimal presenta- 
tions. 




^ closur e ^ 

Figure 5. The closure operation on braids 

The following theorem gives explicit conditions for when two dif- 
ferent braids yield isotopic knots or links. This result was due origi- 
nally to Markov, although a full proof was only published some years 
later by Birman |i5j Theorem 2.3] (see also the paper by Rourke and 
Lambropoulou 1111 ). 

Theorem 1 (Markov LT4J ). Two braids /3i G Bm. and /32 G B„ 
yield closures Pi and P2 which are isotopic as links if and only if 
Pi and p2 are connected by a finite sequence of moves of type Mi 
(conjugation) and M2 (stabilisation), as depicted in Figure^ 
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Figure 6. Markov moves of type Mi (conjugation) and A/2 (stabilisation) 



2 Unknotting 
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Table 1. Braid words for knots with up to eight crossings 1121 



We now attempt to use evolutionary techniques to devise optimal 
unknotting strategies for knots and links. The theorems of Alexan- 
der and Markov enable us to represent knots and links as words in 
the standard generators (t,; of the braid group B„ for some n. The 
crossing-change operation is then simply a matter of taking such a 
word and then replacing a given ai with its inverse or vice- versa. 

Our goal is, given a knot K represented as the closure of a braid 
word w G Bn, to evolve a sequence of certain moves which triv- 
ialises the knot with the smallest number of crossing-changes, thus 
obtaining an upper bound on the unknotting number u{K) of K. The 
allowable moves are those which either leave the isotopy class of the 
corresponding knot unchanged, or which change the sign of a single 
crossing, and may be seen in Table|2] 
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Table 2. Allowable moves on braid words 



In more detail, a move of type cancels a substring of the form 
ata^^, a move of type i?J cancels a substring of the form a~^ai, 
and moves of type introduce corresponding substrings. These are 
the Reidemeister moves of type 2, translated into the context of braid 
words. Similarly, moves of type perform Reidemeister moves of 
type 3, moves of type Mi and A/2 are Markov moves, moves of 
type S represent non-interacting crossings sliding past each other, 
and moves of type U change the sign of a single crossing. 



3 Methods 

In this paper, we use evolutionary techniques to find sequences of 
unknotting primitives which are optimal with respect to two subtly 
different but related problems: 

Problem 1. Given an arbitrary knot K, described as the closure (3 
of some braid word (5 G B„, is there a sequence of moves which 
reduces /3 to the trivial word 1 £ B\? If so, what is the minimal such 
sequence with respect to the number of crossing-change operations 
U^, and in what cases does this yield a sharp upper bound for the 
unknotting number u[K) ofK when compared to known values such 
as those listed in the Knotlnfo database 1121 . 

Problem 2. Given a finite set S of knots, each described as the clo- 
sure /3 of some braid word /3 G Bn, is there a single, universal se- 
quence of moves which, possibly with repeated applications, trivi- 
alises each knot in S ? 



The second of these problems is, broadly speaking, the generali- 
sation of the first to more than one reference knot; equivalently, the 
first problem is the special case of the second where we consider a 
single knot. 



As usual in an evolutionary computation approach, we need to 
define seven things: how population members are represented; the 
parameters; how the population is initialised; how the mutation and 
crossover operators are defined; and how fitness evaluation and se- 
lection happen. 

3.1 Representation 

Each member of the population consists of a list of allowed moves 
drawn from those in Table|2] This list is of variable length. 

3.2 Parameters 

For Problem [T] the population size was 500 and the number of iter- 
ations was 4 X length(_ftr)^, where K is the knot in question and 
length is the number of crossings. For Problem|2]the population size 
was 200 and the number of iterations was 4 x max(length(S'))^, 
where S is the set of knots. These parameters were determined by 
informal experimentation. The program was run three times for each 
knot or set of knots, and the best result is reported. 

3.3 Initialisation 

Each member of the population is initialised by selecting a list of be- 
tween 1 and 15 moves uniformly at random from the moves given in 
Table[2] The moves are selected with replacement; that is, a sequence 
may include more than one move of a given type. 

3.4 Mutation 

There are three different mutation operators, which are selected uni- 
formly at random and applied to individuals with the overall mutation 
probability being 10%. These are: 

1. Select a random move from the list and replace it with a move 
drawn randomly from the moves in Table |2] 

2. Choose a random move from the list and delete it (as long as the 
list of moves contains at least one move). 

3. Choose a random position in the list and insert a randomly chosen 
move from the moves in Table |2] 

3.5 Crossover 

One-point crossover is applied to all individuals as follows: the two 
strings are aligned, a position less than or equal to the length of the 
shortest string is selected at random, and the strings crossed over at 
that point. 

3.6 Fitness evaluation 

The members of the population are evaluated by attempting to unknot 
each of a set of knots, which in the case of Problem [T] will be just a 
single example, and in the case of Problem |2] will consist of more 
than one knot. The execution is carried out as follows. Let M — 
mi, m2, . . . , m„ be the n moves in the list. Let K = Kq be the 
original knot, and Ki, K2, ... be the sequence of knots generated. 

The knot Kq is analysed for the preconditions for mi to be 
carried out, if they are satisfied then the move is applied, so that 
K\ :— mi(A"o); if the preconditions are not satisfied then the knot 
is unchanged (Ki := A'o). The next step is to attempt to apply m2 to 
Ki by seeing if its preconditions are satisfied, and so on. When the 
end of M is reached, the list is begun again from the beginning. 



Each time the knot is changed as a result of applying the move, the 
knot is checked to see if there are any crossings remaining. If so, the 
algorithm terminates, and a positive result returned. 

For Problem[T]we apply the sequence M once only, but for Prob- 
lem|2]we perform repeated applications of M. If the knot hasn't been 
trivialised by 50 applications of A/, we assume that it has become 
stuck in a repeating loop (which, for the 1 701 936 knots with sixteen 
or fewer crossings 1101 . is a valid assumption), terminate the process 
with a negative result, and move on to the next knot (if any) on the 
hst. 

By inspection of Table |2] we see that certain operations (namely, 
those of type R2 and A/2) reduce the length of the braid word under 
investigation, some (types Rs, Mi and U) don't, and the rest (types 
R2 and A/2) increase the length of the braid word. 

It is known (see, for example, the paper by Coward (7)) that there 
exist diagrams for the unknot which can only be reduced to the stan- 
dard (zero-crossing) diagram of the unknot by means of at least one 
move of type i?2. That is, at some point during the reduction process, 
we have to temporarily increase the complexity. In this article how- 
ever, for simplicity, we will restrict ourselves to sequences of moves 
which do not increase the length of the braid word. We intend to 
explore the more general case in later work. 

With regard to Problem[T] we are trying to find a sequence of op- 
erations which unknot a specific knot with the smallest number of 
crossing changes. Slightly less importantly, we want to find the sim- 
plest possible such unknotting sequence. 

For Problem |2] we are trying to find a sequence of operations 
which (perhaps with repeated applications) unknots as many knots 
as possible, as efficiently as possible. 

Let S denote the set of reference braids, on which each sequence is 
being tested. (For Problem[T|this will consist of a single braid word.) 
Let rs{M) denote the number of braids in 5* which are fully reduced 
by (one or more application of) the sequence A/. (In practice, we 
specify an upper threshold of 50 iterations, as described above.) Let 
mins(A/) and maxs(A/) denote, respectively, the minimum and 
maximum number, over all braids in S, of iterations of A/ required 
to reduce a (reducible) braid. Let 1{M) denote the length of the se- 
quence A/. By c(A/) we denote the number of crossing-change (type 
U) operations in A/, and by cg (M) we denote the total number, over 
all braids in S, of successful crossing-change (type U) operations. 
That is, cs{M) gives a measure of the total amount of unknotting 
actually performed by the sequence A/. 

In the case of Problem [1] some of the operations in the string M 
may have no effect on the braid under examination. For example, 
applying an R2 move to a braid which at that stage has no o-^^a^^ 
substrings will leave the braid unchanged, and may thus be safely 
elided from the sequence, resulting in a shorter sequence. Given a 
sequence A/, applied to a braid /3 G /?„, we denote by Zop,(A/) the 
length of the sequence obtained by optimising A/ in this way with 
respect to /3. 

The fitness function should, ideally, seek to minimise the number 
of crossing changes, maximise (at least when working on Problem[2]( 
the number of knots which can be reduced by a given sequence, min- 
imise the length of the sequence, and minimise the number of re- 
peated applications of the sequence required to reduce those braids 
in S which are reducible by the operations under consideration. 

With those criteria in mind, we define the fitness function for Prob- 
lem [T] to be 



/i(A/) = l + 



10000rs(A/) 
/op, (A/) + CS (A/)3-M 



and that for Problem|2]to be 

/2(M) = 1 + 



rs{Mf 



l + maxs(M) + ;(M)' 



Since both /i (M) and /2 (Af) depend only on the set of braids under 
consideration, which doesn't change between generations, we can 
optimise the simulations by caching the fitness values for a given 
string M, rather than recalculating it each time. 

3.7 Selection 

Using an approach similar to the Stochastic Universal Sampling Al- 
gorithm |4|, in each generation, we rank the candidate sequences in 
order of their normalised fitness /(A/): the fitness /(A/) of the se- 
quence M divided by the mean fitness over the whole population. 
The integer part of / (Af) gives the number of copies contributed to 
the next generation, while the fractional part gives the probability 
of an additional copy. So, a sequence M with a normalised fitness 
f{M) = 1.72 contributes one copy to the following generation, plus 
a 72% chance of a second copy. 



3.8 Implementation 

The source code for the implementation is available from the authors 
on request. 

4 Results 

4.1 The single unknotting problem 

TablefTllists braid words for all knots with up to eight crossings, and 
Table [3] lists unknotting sequences for those knots, generated by a 
Perl program implementing Problem [T] The sequences are not nec- 
essarily unique (and in many cases will not be), nor are they guaran- 
teed to be optimal; however we observe that for 21 of the 35 knots 
with eight or fewer crossings, our program has correctly calculated 
the unknotting number. 

Figure [t] shows the braid aio-^^cria^^ (whose closure is iso- 
topic to the figure-eight knot 4i) being reduced by the sequence 
UMiR3R2Ml 




M., 



Figure 7. Reduction of the figure-eight knot 

The knots 3i, 5i and 7i are worth examining a little closer: 
these are the torus knots of type (2, 2n+l), with braid presentation 



on the surface of a standard, unknotted torus so that the strands wind 
p times round the torus in the longitudinal direction, and q times 
in the meridional direction. The torus knot Tp^q of type {p, q) has 
unknotting number «(Tp,q) = |(p—l)(q—l). See Rolfsen JTtI Sec- 
tion 3.C] or Cromwell |8l Section 1.5] for further details on torus 
knots. 
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£ B2 and unknotting number u(K) = n. The unknotting 
sequences obtained for these knots have a very similar form, namely 
{U R2)"~^ M2. More generally, given two positive, coprime integers 
p and q, the torus knot of type {p, q) is the knot which can be drawn 



Table 3. Unknotting sequences M for single knots K, using the braid words 
from Table[T| comparing the number c(M) of crossing changes performed by 
the sequence M, and the crossing number u{K) of the knot K 



4.2 The multiple unknotting problem 

A simulation of Problem |2] again implemented in Perl, obtains uni- 
versal or near-universal unknotting sequences, some examples of 
which are listed in Table|4] Some of the more complex braids (those 
corresponding to knots with minimal crossing number 9 or higher) 
were unreducible by any scheme found by our program, because we 
restricted ourselves to operations which don't increase the length of 
the braid word. As noted earlier (see the paper by Coward 1 7 1 for de- 
tails), sometimes we need to perform a move of type to introduce 
two additional crossings during the reduction scheme. 

5 Conclusions and future work 

In this paper we have seen how an unknotting algorithm can be 
evolved based on a number of primitive moves. There are a number 
of future directions for research in the area of applying evolutionary 
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3i-77 


U H3H20M2AI1 M2 


/: 



1 /I 


14 


3i-82i 


MiUR3M'(R3M2 


5 


35 


35 




SR3R2MI 








3i-82i 


MiUMiR3R2M^ 


7 


35 


35 




R3M1 








3l-949 


UR2R3MIUM2 


6 


74 


84 




R3M1M2M1R3 








3l-949 


R3R2UMlM^ 


10 


73 


84 



Table 4. Universal or near-universal unknotting sequences for multiple 
knots 



algorithms and other machine learning techniques to mathematical 
problems in knot theory and related areas. 

Initially, there are some basic extensions to the work described 
in this paper. For example, rather than focusing on unknotting we 
could use similar techniques to address the related problem of knot 
equivalence. Furthermore, there are similar problems in other areas 
of mathematics (for example, graph theory and group theory) so this 
could be extended to those. 

More interestingly, there is the question of counterexample search. 
There are a number of conjectures in this area where there are some 
measures that could be used to ascertain how close a particular ex- 
ample is to being a counterexample to the conjecture. Using these 
measures, experiments could be done on exploring the space of knots 
to find a counterexample; some preliminary work along these lines 
has been done by Mahrwald 1131 . 

Another approach would be to take a data mining approach to cer- 
tain mathematical problems. For example, we could generate a large 
database of knots and then apply classification techniques to distin- 
guish different classes of knot, or applying clustering techniques to 
group knots according to some metric. An examination of the results 
from this classification might give new insights into the underlying 
structure of the space of knots. A related topic is using genetic pro- 
gramming to evolve invariants, that is, functions that distinguish be- 
tween different knots by processing their diagrams. 
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